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THE PATTERN SPEEDS OF NGC 3031, NGC 2366, AND DDO 154 AS FUNCTIONS OF RADIUS 
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ABSTRACT 

The pattern speeds of NGC 3031, NGC 2366, and DDO 154 are measured using a solution of 
the Tremaine- Weinberg equations derived in a previous paper. Four different data sets of NGC 3031 
produce consistent results despite differences in angular resolution, spectral resolution, and sensitivities 
to structures on different scales. The results for NGC 3031 show that the pattern speed is more similar 
to the material speed than it is to the speed of a rigidly rotating pattern, and that there are no clear 
indications of unique corotation or Lindblad resonances. Unlike NGC 3031, the results for NGC 2366 
and DDO 154 show clear departures from the material speed. The results for NGC 2366 and DDO 
154 also show that the solution method can produce meaningful results that are simple to interpret 
even if there is not a coherent or well-defined pattern in the data. The angular resolution of a data 
set has the greatest affect on the results, especially for determining the radial behavior of the pattern 
speed, and whether there is a single, global pattern speed. 

Key words: galaxies: fundamental parameters - galaxies: individual (NGC 3031, NGC 2366, DDO 154) 
- galaxies: kinematics and dynamics - galaxies: spiral - methods: data analysis - methods: statistical 



1. INTRODUCTION 

The purpose of this paper is to show how different 
properties of a data set affect the results from using 
a method developed in a previous paper (Speights & 
Westpfahl 2011, hereafter Paper I) for solving the pattern 
speed (Q p ) equations of Tremaine & Weinberg (1984, 
hereafter TW84). This is important for choosing data 
in future applications of the solution method, with the 
long-term goal of determining whether spiral patterns 
are rigidly rotating waves, shearing material arms, or 
perhaps some combination of both. 

The pattern speed is traditionally determined from a 
guided interpretation of observations according to the 
density wave theory (Lin & Shu 1964, 1966; Lin et al. 
1969). In density wave theory Q p is a free parameter 
that is measurable by assuming that it is constant (i.e., 
independent of radius) , and that resonant radii are iden- 
tifiable from photometric and kinematic features of the 
galaxy. Resonance is assumed to occur when Q p = S7, 
the material speed (corotation resonance), and when Q p 
— Q ± n/m, where k is the epicycle frequency of a stel- 
lar orbit and m is the mode of the wave corresponding 
to the number of arms (Lindblad resonance). The as- 
sumption that Hp is constant is appealing if one assumes 
that spiral patterns persist for many dynamical (orbital) 
timescales because this avoids the winding dilemma for 
material arms. 

The TW84 equations, however, are independent of the 
detailed assumptions of spiral structure theories, and so- 
lutions can be used to test those assumptions. These 
are integrated forms of the continuity equation that re- 
late the pattern speed to the observable kinematic and 
density properties of a pattern tracer. The derivation 
in TW84 assumes that the disk is flat, that the pattern 
speed is well defined, and that the pattern tracer obeys 
the continuity equation. 

In Paper I the TW84 equations are solved for a pat- 



tern speed that is allowed to vary with radius. The solu- 
tion method is then applied to the Hi spiral pattern of 
the barred, grand-design galaxy NGC 1365. A statistical 
analysis of the results rules out the hypothesis of a con- 
stant Hp for the spiral arms. The pattern speed behaves 
approximately as 1/r, and is very similar to the material 
speed. There are also no clear indications of unique coro- 
tation or Lindblad resonances for the spiral arms. It is 
noted in Paper I that the winding dilemma is avoidable 
for NGC 1365 if the spiral pattern is a transient feature. 

The results for NGC 1365 are inconsistent with tradi- 
tional methods for measuring spiral arm pattern speeds. 
More direct measurements that allow for a radially vary- 
ing flp are therefore needed in order to estimate how 
common shearing patterns are. There is an abundant 
amount of archived data for doing so, and a better un- 
derstanding of how different properties of a data set affect 
the results will help with selecting useful data. 

This paper reports the results from applying the solu- 
tion method of Paper I to the grand-design galaxy NGC 
3031, the irregular galaxy NGC 2366, and the dwarf ir- 
regular galaxy DDO 154. Four different data sets of NGC 
3031 show how angular resolution, spectral resolution, 
and sensitivities to structures on different scales affect 
the results. The pattern speed is found for NGC 2366 
and DDO 154 to show how the results are affected in 
the absence of a coherent or well-organized pattern. The 
results for the irregular galaxies also provide a counterex- 
ample to the results for NGC 3031. 

The rest of this paper is organized as follows. In Sec- 
tion 2 the solution method developed in Paper I is briefly 
reviewed, and a new procedure is introduced for selecting 
models that are an alternative to a constant Q p . In Sec- 
tion 3 the method of analyzing the results is explained. 
In Section 4 the data are described. In Section 5 the re- 
sults are presented and analyzed. In Section 6 the results 
are discussed. In Section 7 is a summary. 
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2. SOLUTION METHOD 

When n p is allowed to vary with radius in the TW84 
equations one of the results is an integral equation, 

f +QC ( d d \ 

J Cl p (r) [x—I(x, y) - y^O, y) ) dx 



+ OC Q 



—I(x,y)v v (x,y)dx, (1) 

that is solvable for simple functional forms of il p . In 
Equation (1), / is the intensity of a pattern tracer and 
v y is the line of sight velocity of the tracer (corrected for 
the systemic velocity, V sys ) times the cosecant of the disk 
inclination, i. The x and y coordinates are measured rel- 
ative to the kinematic minor and major axis respectively. 
For an illustration showing how a galaxy disk is oriented 
in the coordinate system of Equation (1), see Figure (1) 
in Paper I. 

The two forms of fl p that are used in Paper I for solving 
Equation (1) are a polynomial in r = \J x 1 + y 2 , 



i=0 



and a polynomial in 1/r, 



Jh 

«p(r)=E3' 



(2) 



(3) 



i=0 



where that n is an integer and on are unknown coeffi- 
cients that are determined in a solution. Both Equations 
(2) and (3) can represent a constant Sl p when n = 0, and 
higher order terms can be used to determine whether 
the pattern is rigid or shearing. Determining the order 
at which point it is appropriate to truncate the summa- 
tions in Equations (2) and (3) is the primary goal of the 
analysis in Paper I. 

The advantage of using such forms for ft p is that the 
unknown coefficients can come out of the integral on the 
left hand side of Equation (1). Using Equation (2), Equa- 
tion (1) becomes, 

E Q > \x^I{x,y) - y-^I{x,y)} dx 



i=0 



/+oo 
-oo 



d_ 

dy 



I(x,y)v y (x,y) dx, (4) 



or similarly using Equation (3), 

th-L 7i{ x ^ i{x > y) - y -IL i{x > y) } dx 

r+oo q 



dy 



I(x,y)v y (x,y)dx. (5) 



Applying either Equation (4) or (5) independently (i.e., 
at least a beam size apart in the case of radio synthesis 
maps) to different parts of the galaxy transforms the inte- 
gral equation into a matrix equation of the form Ga = d. 
That the calculations are independent is important be- 
cause this is assumed in the methods used for estimating 
and propagating uncertainties, and in the statistics used 



for analyzing the results. The matrix equation is easily 
solved using the normal equations, G T Ga = G T d, for 
the linear least-squares problem. Methods for the anal- 
ysis of the results for linear least-squares problems are 
well established, and this is another advantage of using 
such forms for f2 p . 

The results for NGC 1365 demonstrate that only a few 
terms are needed in the above polynomials in order to de- 
termine the radial behavior of Q p . Including many terms 
can unnecessarily complicate the analysis, and by \n\ = 
3 the matrix G T G typically becomes ill-conditioned for 
numerical inversion. When G T G is ill-conditioned, the 
solution is unstable, and any inferences based on such so- 
lutions may be unreliable. By summing Equations (2) or 
(3) to order n = 3 one can infer whether there is any sign 
of a discontinuity in Q p such as a step function. This is 
useful in Paper I because the bar in NGC 1365 may have 
a separate f2 p than that of the spiral arms (Sellwood & 
Sparkc 1988; Sellwood 1993 and references therein), and 
it needed to be shown whether the pattern speed of the 
bar is detected in the results. 

Since the pattern of NGC 3031 shows no obvious dis- 
continuities or signs of a bar (see Figures 1 through 4), 
and in order to simplify the analysis for all three galaxies, 
only two terms are used in this paper as an alternative 
to a constant fl p model. The alternative models have the 
form, 



(6) 



but without the restriction to integer values of n. Re- 
laxing the restriction to integer values of n allows for 
different amounts of curvature in tt p . The value of n is 
determined by comparing the goodness of fit for solutions 
of Equation (1) using Equation (6) and a range of differ- 
ent values for n. For each solution, the distribution of the 
residuals, normalized by the propagated uncertainties, is 
compared with a standard normal distribution. If one 
assumes that the uncertainties are well known, the dis- 
tribution of the normalized residuals will be a standard 
normal distribution when the true form of Q p is used to 
find a solution. Therefore, the value of n that results 
in a distribution of the normalized residuals that is the 
most similar to a standard normal distribution is used to 
chose an alternative model for subsequent analysis. 

A Kolmogorov - Smirnov (KS) test is used to quantify 
how similar the distribution of the normalized residu- 
als is to a standard normal distribution. The KS test 
consists of finding the largest distance between two cu- 
mulative distribution functions (cdf), D = Max|cdf(i?) 
— cdf (SN) | , where that in this case R is the distribution 
of the normalized residuals and SN is a standard nor- 
mal distribution. Then, the probability, Pks, is found 
for obtaining a value of D should the null hypothesis be 
true that R is a standard normal distribution (Kendall 
& Stewart 1973, Chapter 30; Lupton 1993, Chapter 19). 
If Pks is sufficiently small, then there is convincing evi- 
dence that R ^ SN, and one can infer that the model in 
question is not the true model. The following convention 
is adopted for interpreting values of Pks- Pks < 1% is 
convincing evidence that R ^ SN; Pks ~ 5% is sug- 
gestive of R ^ SN, but inconclusive; and larger values 
cannot rule out the null hypothesis that R is a standard 
normal distribution. In this paper a lower limit of 0.01% 
is chosen as a cutoff for reporting small probabilities. As 
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a compliment to the analysis explained in the next sec- 
tion, the KS test is also applied to the results from using 
a constant f2 p model. 

For the purpose of determining n, a larger value of 
Pks is considered an indication that the distribution of 
the normalized residuals is more similar to a standard 
normal distribution. In Section 5 it is shown that two 
peaks consistently occur in plots of Pks vs n i one each 
for negative and positive values of n. The results are 
analyzed for each case, where the model using a negative 
value of n = n\ is referred to as alternative model 1 
(AMI), and the model using a positive value of n = n-2 
is referred to as alternative model 2 (AM2). 

3. ANALYSIS 

The goal of the analysis is to determine whether the 
pattern is rigid or shearing. This is based on a compar- 
ison of reduced chi-square (xt f° r v degrees of freedom) 
for a solution using a constant fl p model with the values 
obtained for AMI and AM2. Unless otherwise noted in 
Section 5, the validity of using xt as an indicator of the 
goodness of fit is based on the results of a Lilliefors test 
for quantifying how similar R is to a normal distribu- 
tion. The Lilliefors test follows the same procedure as 
the KS test, with the difference that D = Max|cdf(i?) 
— cdf (N) | , where that N is a normal distribution with a 
mean and a standard deviation that is the same as that 
for R (Mandansky 1988, Chapter 1). If the probability, 
Pl, of obtaining D is sufficiently small, then R may not 
have a normal distribution, and the value of xt m ay not 
be a reliable indicator of the goodness of fit. The same 
convention that is adopted for interpreting values of Pks 
is adopted for interpreting values of Pl ■ 

If a solution using AMI or AM2 has a value of \t that 
is an improvement to the value of \t f° r the constant Q p 
model, then Fisher's F test of an additional term is used 
to determine if the improvement is statistically signifi- 
cant. The F test is a natural choice for quantifying the 
significance of the improvement in xt because the con- 
stant fi p model is a reduced version of AMI and AM2. 
The F test consists of calculating F — A\ 2 /xt an d find- 
ing the probability, Pp, for obtaining such a value of F 
should the null hypothesis be true that the extra term, 
«i, is zero (Bevington & Robinson 2002, Chapter 10). 
If Pp is sufficiently small, then there is convincing evi- 
dence that a± ^ 0, and one can infer that the pattern is 
shearing. If such evidence is found then the result for the 
constant Q p model is the mean of f2 p . The following con- 
vention is adopted for interpreting values of Pp: Pp < 
1% is convincing evidence that a± ^ 0; Pp ~ 5% is sug- 
gestive of a% =fi 0, but inconclusive; and larger values are 
an indication that the improvement is not statistically 
significant. 

The uncertainties for calculating xZ' confidence inter- 
vals for the coefficients, and confidence bands for O p , are 
estimated by propagating the uncertainties per pixel in 
the maps of I and v y , as well as the uncertainties in x, 
y, i, and V S y S through the calculations of both G and d. 
The method of estimating the uncertainties in the maps 
of / and v y are adopted from the reasoning given in Pa- 
per I. This takes into account the rms noise, the channel 
width, and the number of channels a typical detection 
occurs. The total uncertainty for the jth independent 



calculation of Equation (1) is cr| — a\. + Y17=o a G i a l 
(Bevington & Robinson 2002, Chapter 6). The confi- 
dence intervals and confidence bands are not estimated 
using the covariance matrix because of the uncertainties 
in G. Instead, these are estimated using Monte Carlo 
methods. 

An incorrect value for the position angle of the disk 
(PA) can potentially cause large errors in solutions of 
the TW84 equations (Debattista 2003; Debattista & 
Williams 2004). The PA does not explicitly appear in 
Equation (1), and consequently the uncertainty in the 
PA is not accounted for in the calculation of the propa- 
gated uncertainties. How an incorrect PA could affect the 
conclusions in this paper is investigated in Section 5 by 
changing the PA and applying the solution method and 
analysis. This is done for a range of -2er PA ^ APA ^ 2cr PA 
where that cr PA is the la uncertainty for the adopted PA. 
In Paper I it is explained that this method is preferred 
to Monte Carlo methods because Monte Carlo methods 
typically produce non-gaussian results for the distribu- 
tions of the coefficients. 

In order to minimize the uncertainties in the analysis, 
instrumental units of [fl p ] — km s -1 arcsec -1 are used for 
solving Equation (1) and analyzing the results. This is 
preferred to units commonly used for galactic disks, [Cl p ] 
= km s _1 kpc -1 , because the conversion from arcseconds 
to kpc adds a large amount of uncertainty to x and y 
owing to the uncertainty in the distance to the galaxy. 
Furthermore, units commonly used for galactic disks are 
not necessary for the purpose of this paper. 

4. DATA SETS 

Radio synthesis intensity and velocity maps of H I are 
used for applying the solution method. For a discussion 
of the advantages of using H I as a pattern tracer the in- 
terested reader is referred to Paper I. Although TW84 
use a form of the continuity equation that assumes mass 
conservation, the TW84 equations are derivable form the 
general form of the continuity equation that includes a 
source function for the creation and destruction of mate- 
rial (Westpfahl 1998). If the source function is constant 
within dynamical timescales, then it does not affect the 
calculation of the pattern speed (Paper I) . For the pur- 
pose of this paper it is assumed that the source function 
is approximately constant within dynamical timescales. 

The disk parameters that are adopted for measuring 
fi p are provided in Table 1. The values for the kinematic 
centers are adopted from Trachternach et al. (2008). The 
values for i, PA, and V S ys are adopted from de Blok et 
al. (2008) for NGC 3031, and from Oh et al. (2011) for 
NGC 2366 and DDO 154. 

Figures 1 through 6 show the maps of / and v y for 
the four different data sets of NGC 3031, and the data 
sets of NGC 2366 and DD0154, respectively. The origi- 
nal maps were modified into their current form using the 
Astronomical Image Processing System developed and 
maintained by the National Radio Astronomy Observa- 
tory (NRAO). These are rotated and shifted from their 
original orientation so that their kinematic major axes 

1 The National Radio Astronomy Observatory is operated by 
Associated Universities, Inc., under cooperative agreement with 
the National Science Foundation. 
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TABLE 1 

Adopted Disk Parameters 



Galaxy 


Kinematic Center 




PA 

(°) 






H. A. (J2000) 

m s \ 


Dec. (J 2000) 
(° ' ") 


i 

(°) 


Vsys 
(km s ) 


NGC 3031 


09 55 33.5 ± 0.6 


+69 03 52.0 ± 3.9 


59 ± 2.6 


330 ± 1.3 


-39 ± 2.8 


NGC 2366 


07 28 53.9 ± 0.7 


+69 12 37.4 ± 7.8 


63 ± 1.6 


39 ± 4.0 


104 ± 2.0 


DDO 154 


12 54 05.9 ± 0.2 


+27 09 09.9 ± 3.4 


66 ± 2.5 


229 ± 7.8 


375 ± 1.5 
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Fig. 1. — Data for NGC 3031 from AW96a. Panel (a) shows is an optical image at the same scale for comparison. Panels (b) and (c) 
show the maps of / and v y , respectively. Directions of north (N) and east (E) on the sky are indicated in panel (a). The receding half of 
the galaxy is to the right. The peak flux is 134 Jy beam -1 m s . The velocities are binned in increments of 67 km s _1 . 
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Fig. 2. — Data for NGC 3031 from AW96b. The figure is formatted in the same way as Figure 1. The peak flux is 2,880 Jy beam -1 m 
s — 1 . The velocities are binned in increments of 67 km s — 1 . 



are along the x axis, and the origins of the x-y coor- 
dinates are the kinematic centers that are adopted for 
each galaxy. Included in the figures are optical images 
from the Digitized Sky Survey. These are based on the 
photographic data of the National Geographic Society - 
Palomar Observatory Sky Survejd (NGS-POSS), and are 



provided for comparison with the Hi maps. The optical 
images are from the 48 inch Samuel Oschin Telescope 
on Palomar Mountain and a 103aE (6450 A) emulsion. 
Table 2 provides a summary of the properties of each 
data set. The data are described in greater detail in the 
following subsections. 



2 The NGS-POSS was funded by a grant from the National Ge- 
ographic Society to the California Institute of Technology. The 
plates were processed into the present compressed digital form 



with their permission. The Digitized Sky Survey was produced 
at the Space Telescope Science Institute under US Government 
grant NAG W-2166. 
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TABLE 2 

Properties of the Data 



Galaxy 


Notes 


Channel Seperation 


Beam FWHM 


rms 






(km s" 1 ) 


(") 


(mjy beam -1 ) 


NGC 3031 


AW96a 


2.6 


9.0 x 9.0 


0.6 


NGC 3031 


AW96b 


2.6 


60.0 x 60.0 


1.23 


NGC 3031 


H84 


10.3 


9.0 x 9.0 


2.5 


NGC 3031 


YP99 


6.6 


135.3 x 124.8 


6.8 


NGC 2366 




2.6 


13.1 x 11.9 


0.5 


DDO 154 




2.6 


14.1 x 12.6 


0.5 
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Fig. 3. — Data for NGC 3031 from H84. The figure is formatted in the same way as Figure 1. The peak flux is 497 Jy beam 1 
The velocities are binned in increments of 67 km s _1 . 
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Fig. 4. — Data for NGC 3031 from YP99. The figure is formatted in the same way as Figure 1. The peak flux is 2.291 Jy beam -1 m 
s~ 1 . The velocities are binned in increments of 67 km s — 1 . 



4.1. Data Sets for NGC 3031 

The four different data sets of NGC 3031 are from three 
different observations (Adler & Westpfahl 1996, hereafter 
AW96; Hine 1984, hereafter H84; Yun & Purton 1999, 
private communication, hereafter YP99). They have dif- 
ferent combinations of angular resolution, spectral reso- 
lution, and sensitivities to structures on different scales. 

Two of the four data sets, hereafter referred to as 
AW96a and AW96b, are made from the observations 



of AW96. They used the Very Large Array (VLA) of 
the NRAO in 1992 and 1993 with all 27 antennas, 128 
channels, one intermediate frequency (IF), and on-line 
Hanning smoothing. Data from two pointings of the 
VLA in the B, C, and D configurations with a chan- 
nel width and separation of 12.2 kHz (2.6 km s _1 ) are 
combined to make a mosaic. The maps published by 
AW96 do not have the continuum emission of the point 
sources subtracted. For the AW96a data set this emis- 
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Fig. 5. — Data for NGC 2366. The figure is formatted in the same way as Figure 1. The peak flux is 1,020 Jy beam -1 m s -1 . The 
velocities are binned in increments of 15 km s . Note the non-circular velocities for the lower left quadrant of the HI disk. 




Fig. 6. — Data for DDO 154. The figure is formatted in the same way as Figure 1. The peak flux is 455 Jy beam 1 m s . The velocities 
are binned in increments of 14 km s . 



sion is subtracted from a calibrated and mapped data 
cube provided by Adler (2010, private communication). 
The continuum subtracted data cube was used to make 
new intensity and velocity maps. For this data set the 
rms signal in a line-free channel is 0.6 mJy beam -1 , and 
the FWHM of the synthesized beam is a 9.0" circle. The 
AW96b data set is from the same calibrated u-v data, 
also provided by Adler (2010, private communication). 
The u-v data for this data set is tapered using a gaus- 
sian weighting function so that the synthesized beam is a 
60.0" circle. For the AW96b data set the rms signal in a 
line-free channel is 1.23 mJy beam -1 . The details of the 
observing and data reduction are described by AW96. 

The data of H84 are made from observations using the 
VLA in 1982. Owing to limitations of the telescope and 
correlator at that time only 18 antennas and 32 spec- 
tral channels could be used with one IF. Data from two 
pointings of the VLA in the B and D configurations with 
a channel separation of 48.8 kHz (10.3 km s -1 ) are com- 
bined by H84 to make a mosaic. The rms signal in a 



line- free channel is 2.5 mJy beam -1 . The FWHM of the 
synthesized beam is a 9.0" circle. This data set includes 
the same range of baselines as does the AW96 data sets, 
but with relatively fewer intermediate ones due to the 
exclusion of the C configuration. The details of the ob- 
serving and data reduction are described by H84. 

The data of YP99 are made from observations us- 
ing the synthesis telescope of the Dominion Radio As- 
trophysical ObservatorjQ in 1993. The calibrated and 
mapped data cube is kindly made available in advance 
of publication. The channel separation is 15.63 kHz (3.30 
km s -1 ), the channel width is 24.99 kHz (5.28 km s -1 ). 
Pairs of channels are averaged to improve the signal-to- 
noise ratio, giving final channel separations of 31.25 kHz 
or 6.60 km s . The rms signal in a line-free channel 
of this width is 6.8 mJy beam -1 . The FWHM of the 
synthesized beam is 135.3" x 124.8" with a beam PA of 
-4.4°. Regions of statistically significant signal in this 

3 The Dominion Radio Astrophysical Observatory is operated as 
a national facility by the National Research Council of Canada. 
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cube are identified by first convolving each channel with 
a circular Gaussian of diameter 200". Regions in the 
convolved channel with signal of absolute value less than 
2.5 mjy beam -1 are judged to contain no significant de- 
tection, and these regions are blanked in the original, 
unconvolved channel. The blanked channels are used to 
make the data maps. 

The different data sets emphasize different attributes 
of NGC 3031. The YP99 data set excels at showing 
the faint outer structure, but its large synthesized beam 
cannot show fine detail. The AW96a data set excels at 
showing fine detail, but AW96 show that it is missing 
flux, particularly on large scales. The AW96b data set 
does not show the fine detail of the AW96a data set due 
to the u-v tapering. The H84 data set is in some ways 
a compromise between the AW96a and YP99 data sets. 
It does not have the sensitivity to the small structures 
of AW96a, but it is missing little or no flux from NGC 
3031 itself. More fine detail is shown in Figure 1 for 
the data set from AW96a than in Figure 3 for the data 
set from H84 despite having the same sized synthesized 
beam. This is because of the much higher sensitivity of 
the AW96a data set.. 

4.2. Data Sets for NGC 2366 and DDO 154 

The data sets for NGC 2366 and DDO 154 are from 
the THINGS Hi survey (Walter et al. 2008). Both data 
sets were made from B, C, and D observations using the 
VLA. The observations of NGC 2366 are made in 2003 
and 2004, and of DDO 154 in 2004 and 2005. Both data 
sets have a channel width of 2.6 km s" 1 and an rms signal 
in a line- free channel of 0.5 mjy beam" 1 . The FWHM of 
the synthesized beam is 13.1" x 11.9" with a beam PA of 
-1.6°' for NGC 2366, and 14.1" x 12.6" with a beam PA 
of -34.0° for DDO 154. The details of the observing and 
data reduction are described by Walter et al. (2008). 

5. RESULTS 

The calculations for solving Equation (1) and analyz- 
ing the results are carried out in MATLAE0- Differentia- 
tion is performed using the Savitzky-Golay method (Sav- 
itzky & Golay 1964). Integration is performed by sum- 
ming over a row of pixels in a map of the integrand. The 
KS test and the Lilliefors test are performed using MAT- 
LAB's Statistics Toolbox. The rest of the beginning of 
this section explains how the results are presented. The 
details of the results and their analysis are discussed in 
the subsections that follow. 

Plots of Pks vs n are shown in Figure 7. Note that 
two peaks consistently occur in the plots. It is assumed 
that rip does not have a large amount of curvature so 
that ni is not <C 1 and is not 3> 1. Out of a range 
of -6 ^ n ^ 6 for which solutions are found, the peaks 
occur within -2 ^ n ^ 2, so only this range is shown in 
the figure. The values of Pks are not shown in the figure 
for when n = because the determinant of G T G is zero 
for a model of fl p = cxq + otx. The values of n, D, and 
Pks for the peaks in Figure 7 are summarized in Table 
3. Values of D and Pks for the constant fl p model are 
also shown in Table 3. The details of the results for the 

4 The Mathworks, Version 7.10.0.499 (R2010a), 
http: / /www. mathworks. com 



KS test are discussed in paragraph 2 of Sections 5.1, 5.2, 
and 5.3. 

The results for the coefficients are shown in Table 4. 
The uncertainties for the coefficients are reported as the 
half-widths (HW) of their 95% confidence intervals. The 
coefficients are used to calculate plots of Q p for compari- 
son with f2 and possible locations for Lindblad resonance. 
The plots of ftp calculated from the coefficients are dis- 
cussed in paragraph 3 of Sections 5.1, 5.2, and 5.3. 

The results for the F test are shown in Table 5. An 
asterisk next to a value of xt indicates that the normal- 
ized residuals either failed a Liliefors test for normality, 
or that the results of the test are inconclusive. The la 
uncertainties for xt are estimated from the variance, 2v, 
of a x 2 distribution. These show how similar the values 
of xl are for a given pair of AMI and AM2. The details 
of the results for the F test are discussed in paragraph 4 
of Sections 5.1, 5.2, and 5.3. 

5.1. Results for NGC 3031 

Results are presented for the entire spiral pattern, 
which consists of the region \y\ ^ 600" in the Hi maps. 
This provides 42 independent calculations of Equation 
(1) for the AW96a and H84 data sets, 10 independent 
calculations for the AW96b data set, and 4 independent 
calculations for the YP99 data set. As a check, solutions 
are also found for three other regions: \y\ ^ 300", —600" 
^ y sC 0", and 0" sC y < 600". The results for the three 
other regions are consistent with those for the entire spi- 
ral pattern, so they are not presented separately. 

In Figure 7 there are prominent peaks in the values 
of Pks for all four data sets. The values of Pks are 
generally quite large at their peaks, with the exception 
of the peaks at ri2 for the AW96b and YP99 data sets. 
The AW96a, AW96b, and H84 data sets show small P K s 
near n = 0, but this is not the case for the YP99 data set. 
For all four data sets the values of Pks for the constant 
fl p model are small enough to rule out the null hypothe- 
sis that the normalized residuals have a standard normal 
distribution. This is convincing evidence that the con- 
stant rip model is not the true model. This conclusion, 
however, is only reliable for the AW96a, AW96b, and 
H84 data sets because the KS test becomes unreliable 
when there are only a few calculations of the normalized 
residuals. This may explain the discrepancy between the 
values of Pks near n = in Figure 7, and the small 
value of Pks f° r the constant fi p model in Table 3 for 
the YP99 data set. 

Plots of fl p calculated from the coefficients for NGC 
3031 are shown in Figure 8. Included in the plots is 
the material speed, il, and possible locations for 2-arm 
Lindblad resonance. The material speed is calculated 
from the rotation curve by de Blok et al. (2008). The 
figure shows consistent plots among the four data sets 
for each of the different forms of fl p . The plots for AMI 
show the same radial behavior as fl, with smaller val- 
ues of ftp in the inner half of the disk for the AW96b 
data set, and outer half of the disk for the YP99 data 
set. The solutions using n\ more closely resemble the 
concave up behavior of fl in the inner half of the disk 
whereas solutions using ri2 more closely resemble the lin- 
ear behavior of f2 in the outer half of the disk. The plots 
of AMI and AM2 for the AW96a, AW96b, and H84 data 
sets do not show a unique point of corotation resonance. 
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TABLE 3 



Results for the KS Test 


(~Xp\ lavv 

\J CXI CXjV Y 


Notes 


Constant fl p 




AMI 






AM2 








D 




m 


D 


f KS 


ri2 


D 








x 10" 1 


(%) 




x 10" 1 


(%) 




x 10" 1 


(%) 


NGC 3031 


AW96a 


3.64 


< 0.01 


-0.89 


1.18 


56.4 


0.80 


1.46 


30.5 


NGC 3031 


AW96b 


6.00 


0.06 


-0.80 


1.52 


94.9 


0.36 


3.44 


14.7 


NGC 3031 


H84 


4.15 


< 0.01 


-0.61 


0.94 


82.1 


0.84 


0.84 


90.6 


NGC 3031 


YP99 


7.50 


0.78 


-1.42 


2.47 


90.9 


0.31 


4.99 


18.9 


NGC 2366 


A 


3.33 


0.45 


-0.11 


1.09 


88.4 


0.12 


1.01 


93.2 


NGC 2366 


B 


6.14 


< 0.01 


-0.06 


1.58 


82.3 


0.08 


1.55 


84.2 


NGC 2366 


C 


4.02 


2.09 


-0.20 


1.38 


93.5 


0.25 


1.45 


91.2 


NGC 2366 


V 


3.32 


8.89 


-0.31 


2.09 


55.2 


0.41 


1.52 


88.0 


DDO 154 




1.79 


37.9 


-1.01 


1.62 


50.4 


1.82 


2.11 


20.6 




Fig. 7. — Plots of the results from the KS test. The solid lines are values of Pks- The dashed lines at Pks = 5% are provided for 
reference. Panels (a), (b), (c), and (d) show the results for the AW96a, AW96b, H84, and YP99 data sets of NGC 3031, respectively. 
Panels (e), (f), (g), and (h) show the results for regions A, B, C, and T> of NCG 2366, respectively (See Section 5.2 for an explanation of 
these regions). Panel (i) shows the results for DDO 154. 



For the same plots the existence of an outer 2-arm Lind- 
blad resonance is inconclusive because the radius where 
fi p intersects possible locations for such a resonance also 
show corotation. For both the AW96a and H84 data 
sets, the plots of AMI and AM2 do not show an inner 
2-arm Lindblad resonance. The existence of an inner 
2-arm Lindblad resonance is not ruled out by AM2 for 



the AW96b data set due to the slightly larger confidence 
bands. The existence or absence of these resonances in 
plots of AMI and AM2 are not as well constrained for 
the YP99 data set due to the larger confidence bands. 

For all four data sets, and both AMI and AM2, the val- 
ues of Pp are small enough to rule out the null hypothesis 
that a\ = 0. This is convincing evidence that the pat- 
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TABLE 4 

Results for the Coefficients 



Galaxy Notes Constant Cl p AMI AM2 



a ± HW a ± HW 51 ± HWi a ± HW <*i ± HWi 

(km s _1 arcscc -1 ) (km s _1 arcscc -1 ) (km s _1 ) (km s _1 arcscc -1 ) (km s _1 arcsec -2 ) 



NGC 3031 


AW96a 


(4.02 ± 0.51) 


X 


io- 


l 


(-1.20 ± 0.53) 


X 


io- 


i 


(1.50 ± 


0.15) 


X 


10 2 


0.95 ± 0.09 


(-0.32 ± 0.54) 


X 


10" 


-2 


NGC 3031 


AW96b 


(2.12 ± 1.20) 


X 


io- 


i 


(-1.04 ± 0.64) 


X 


io- 


i 


(0.73 ± 


0.13) 


X 


10 2 


1.13 ± 0.24 


(-7.19 ± 1.84) 


X 


10" 


-2 


NGC 3031 


H84 


(3.84 ± 0.60) 


X 


io- 


l 


(-4.57 ± 1.70) 


X 


io- 


l 


(0.49 ± 


0.10) 


X 


10 2 


1.07 ± 0.16 


(-0.25 ± 0.06) 


X 


io- 


-2 


NGC 3031 


YP99 


(3.11 ± 3.11) 


X 


io- 


l 


(-0.55 ± 4.20) 


X 


io- 


i 


(25.21 ± 26.49) 


X 


10 2 


1.35 ± 0.84 


(-14.7 ± 11.8) 


X 


io- 


-2 


NGC 2366 


A 


(2.32 ± 0.73) 


X 


io- 


l 


-1.48 ± 0.82 








3.00 ± 


1.44 






1.71 ± 0.72 


(-7.76 ± 3.76) 


X 


io- 


-1 


NGC 2366 


B 


(2.66 ± 1.06) 


X 


io- 


l 


-2.63 ± 2.96 








3.93 ± 


4.02 






2.51 ± 1.84 


(-15.3 ± 12.5) 


X 


10" 


-1 


NGC 2366 


C 


(2.13 ± 1.00) 


X 


io- 


l 


-1.09 ± 0.69 








3.60 ± 


1.89 






1.21 ± 0.53 


(-2.66 ± 1.42) 


X 


io- 


-1 


NGC 2366 


V 


(2.63 ± 1.24) 


X 


io- 


l 


-0.21 ± 0.44 








2.29 ± 


2.09 






0.61 ± 0.35 


(-0.39 ± 0.38) 


X 


io- 


-1 


DDO 154 




(1.76 ± 0.33) 


X 


io- 


l 


(-1.22 ± 0.29) 


X 


io- 


l 


(1.26 ± 


0.46) 


X 


10 1 


(2.15 ± 0.50) x 10- 1 


(-1.33 ± 1.36) 


X 


10" 


-6 



TABLE 5 

Results for the F Test 



Galaxy 


Notes 


Constant Q p 




AMI 






AM2 








xi 


± a 


xi ± ° 


F 


Pf 

(%) 


xi ±° 


F 


Pf 

(%) 


NGC 3031 


AW96a 


25.8 


± 0.22 


1.16 ± 0.22 


873 


< 0.01 


1.06 ± 0.22 


955 


< 0.01 


NGC 3031 


AW96b 


125 


± 0.47 


4.28 ± 0.50 


255 


< 0.01 


2.12 ± 0.50 


523 


< 0.01 


NGC 3031 


H84 


24.0 


± 0.22 


1.29 ± 0.22 


874 


< 0.01 


1.09 ± 0.22 


1036 


< 0.01 


NGC 3031 


YP99 


73.0 


± 0.82 


3.40 ± 1.00 


52.8 


1.84 


2.19 ± 1.00 


97.9 


1.01 


NGC 2366 


A 


9.54 


± 0.28 


0.90 ± 0.29 


241 


< 0.01 


0.87 ± 0.29 


250 


< 0.01 


NGC 2366 


B 


6.97 


± 0.39* 


0.65 ± 0.41 


127 


< 0.01 


0.83 ± 0.41 


97.6 


< 0.01 


NGC 2366 


C 


13.8 


± 0.41 


1.02 ± 0.43 


151 


< 0.01 


1.09 ± 0.43 


141 


< 0.01 


NGC 2366 


V 


5.55 


± 0.41 


3.52 ± 0.43* 


7.94 


1.67 


3.07 ± 0.43 


10.7 


0.75 


DDO 154 




0.83 


± 0.29 


0.65 ± 0.30 


6.71 


1.67 


0.39 ± 0.30 s 


26.2 


< 0.01 



tern is shearing. Note that this conclusion is consistent 
with the results of the K S test applied to the constant 
Vt p model. The results for the constant fi p models are 
therefore estimates of the mean of il p . For each data set 
it is not clear whether AMI or AM2 provides the bet- 
ter solution, and thus the better estimate of il p . Their 
values of xi are quite similar, especially given their es- 
timated la uncertainties. Despite this ambiguity, it is 
possible to infer from the the above conclusions and the 
plots of AMI and AM2 that there are no clear signs of 
2-arm Lindblad resonances or a unique location of coro- 
tation resonance for the spiral pattern of NGC 3031. It is 
also possible to infer that the speed of the spiral pattern 
more closely resembles the speed of the material than the 
speed of a rigidly rotating pattern. 

Figure 9 shows how an incorrect PA could affect the 
conclusions for the NGC 3031 data sets. In the figure 
(and similar figures for NGC 2366 and DDO 154) are 
values of n, D, xii and the mean of Q p that result 
from using a range of PA within ± 2a PA of the assumed 
PA. The deviations for the mean of fl p are shown because 
they are indicative of the amount of deviation in AMI 
and AM2. The values of n and xi shown in the figure for 
AMI and AM2 do not vary as much for the AW96a and 
H84 data sets as they do for the AW96b and YP99 data 
sets. For all four data sets the values of D are too small to 
rule out the null hypothesis that the normalized residuals 
for AMI and AM2 have a standard normal distribution, 
except for the H84 data set at APA ~ -1.5°. This means 
that the AMI and AM2 for the H84 data sets may not be 
the true model if the assumed PA is incorrect by ~ 1.5°. 
Determining whether the pattern is shearing, however, 
does not require that AMI or AM2 are the true model. 
A APA as large as ± 2er PA would not change the results 



of the F test for the AW96a, AW96b and H84 data sets. 
Although some of the AMI F values around APA < 1° 
for AW96b appear close to the value at which Pp = 5% 
in the figure, none of the F values are actually small 
enough that Pp ^ 1%. There are 6 points for the YP99 
data set where that the F test is inconclusive. These 
occur at APA = -0.94° for AMI and -0.73° < APA «C 
-0.31° for AM2. The only point in Figure 9 where that 
the F test cannot rule out the null hypothesis that a\ 
= occurs at APA = -0.83° for the YP99 data set and 
AMI. 

Two points are worth emphasizing regarding how an 
incorrect PA could affect the conclusions for NGC 3031: 

1. ) All of the values of F shown for three out of four 

data sets arc large enough to conclusively rule out 
the null hypothesis that ct\ = 0. For all four data 
sets, 98.5% of the F values are large enough to 
conclusively rule out the null hypothesis. 

2. ) Most of the mean fi p results (87.7%) are within the 

95% confidence bands for the mean of Q p obtained 
using the adopted PA. The most deviant points 
occur for the AW96b data set, and these are larger 
than this data set's mean value for f2 p . Such an 
increase in the fl p results for the AW96b data set 
would only bring the results for this data set into 
better agreement with the results for the other data 
sets. 

Therefore, the affect of an incorrect PA on the conclu- 
sions for NGC 3031 are insignificant. The affect of an 
incorrect PA on the conclusions for NGC 1365 were also 
found to be insignificant in Paper I. The robustness of 
the conclusions to an incorrect PA, in both cases, is at- 
tributable to the small uncertainties in the measured val- 
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Fig. 8. — Plots of f2 p and possible locations for resonance in NGC 3031. The solid lines (colored red in the online version) and dashed 
lines show Q p and the 95% confidence bands, respectively. The open circles (filled in green in the online version) show the material speed, 
Q. The dash-dot lines (colored blue in the online version) show possible locations of 2-arm Lindblad resonance. From the top row to the 
bottom, shown are the results for the AW96a, AW96b, H84, and YP99 data sets. The left, center, and right columns of panels show the 
results for a constant f2 p , AMI, and AM2, respectively. Note how similar the plots for the AMI and AM2 models are to Q. 



ues of the PA, and the high quality of the Hi data. 

5.2. Results for NGC 2366 

Unlike NGC 3031, different regions of NGC 2366 pro- 
duce different results. Results are presented for four dif- 
ferent regions: \y\ 400", \y\ ^ 200", -400" < y ^ 0", 
and 0" ^ y ^ 400". These are hereafter referred to as 
regions A, B, C, and V, respectively. They provide 26, 
14, 13, and 13 independent calculations of Equation (1), 
respectively. 

In Figure 7 there are prominent peaks in the values of 



Pks for regions A, B, and C. The peaks for region V are 
not as prominent. Similar to the plots for NGC 3031, the 
plots for all four regions of NGC 2366 show small Pks 
near n = 0. For regions A, B, and C, the values of Pks 
for the constant Q p model are small enough to rule out 
the null hypothesis that the normalized residuals have 
a standard normal distribution. This is convincing evi- 
dence that the constant Q p model is not the true model 
for these regions. For region £>, Pks for the constant f2 p 
model is too large to rule out the null hypothesis that 
the normalized residuals have a standard normal distri- 
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Fig. 9. — Plots showing how an incorrect PA could affect the findings for the NGC 3031 data sets. From left to right are sets of plots 
for the AW96a, AW96b, H84, and YP99 data sets. From top to bottom, each set of plots shows values of n, D, and F where that the 
values for AMI and AM2 are shown as open squares (colored red in the online version) and filled in diamonds (colored blue in the online 
version), respectively. The bottom plot in each set shows the mean of Q p as open circles (colored blue in the online version). The dashed 
lines (colored green in the online version) included in the plots of D and F show the value where Pp> and Pp are equal to 5%. Note that 
smaller values of D correspond to larger values of Pjj, and vice versa for F and Pp. The solid lines (colored red in the online version) 
and dashed lines included in the plots of the mean of f2 p are the values obtained using the assumed PA and their 95% confidence bands, 
respectively. 



bution. Nothing can be inferred from the KS test about 
the constant f2 p model in region T>. 

Plots of Sip calculated from the coefficients for NGC 
2366 are shown in Figure 10. The material speed in the 
figure is calculated from the rotation curve by Oh et al. 
(2008; hereafter O08). Possible locations for 2-arm lind- 
blad resonances are included in the plots (and also for 
DDO 154) to create the same sense of scale as the plots 
for NGC 3031. All of the regions show consistent plots 
for the constant Q, p model. All of the regions also show 
corotation in the inner half of the disk for AMI and AM2. 
Note that for each region, their plots for AMI and AM2 
are similar although their values of n have opposite sign. 
There are clear differences among the different regions 
for the plots in the outer half of the disk. These are the 



most prominent for regions C and T>. In region C the 
plots for AMI and AM2 approach fl p = at ~ 400", 
which is ~ 0.16 km s _1 arcsec -1 less than Q. In region 
T> the results for AMI and AM2 are similar to fi for the 
entire range in r shown. 

For all four regions, and both AMI and AM2, the val- 
ues of Pp are small enough to rule out the null hypothesis 
that ot\ — 0. For regions A, B, and C, the results of the 
F test are consistent with the results of the KS test ap- 
plied to the constant O p model. The results of the F test 
for region B, however, are unreliable because the distri- 
bution of the normalized residuals for the constant £l p 
model failed a Lilliefors test for normality (Pl =1.79%). 
The F test result for AMI in region T> may also be unre- 
liable because the results of a Lilliefors test applied to its 
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Fig. 10. — Plots of Q p and possible locations for resonance in NGC 2366. The figure is formatted in the same way as Figure 8. From the 
top row to the bottom, shown are the results for regions A, B, C, and T>. The left, center, and right columns of panels show the results for 
a constant f2 p , AMI, and AM2, respectively. Note the differences for the results shown in panels (h) and (fc), and likewise for panels (i) 
and (I). 



normalized residuals is inconclusive (Pl — 5.35%). From 
the results for the F test, the results of the KS test for 
the constant O p model, and the plots of Q p , it is possible 
to infer that the clumps of H I in NGC 2366 are corotat- 
ing with the material in the inner half of the disk, and 
that there are clumps of H I in the outer half of the disk 
in region C that are rotating more slowly than the rest 
of the material. 

The results for NGC 2366 show that the solution 
method can detect different values of fl p in different re- 
gions of a galaxy. This demonstrates the usefulness of 



the solution method for determining if there is a single, 
global pattern speed. The difference for regions C and 
T> may be attributed to the non-circular motions in the 
velocity field that are present in region C. These are eas- 
ily identified in Figure 2 by the kinks in the contours of 

constant v y around y ~ 200 - 400 . The non-circular 
motions were removed from the velocity field when O08 
derived the rotation curve used in this paper for calcu- 
lating O. The difference between the rotation curve of 
O08, and one obtained without the non-circular motions 
removed (see Figure 13 of O08), is not large enough to 
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Fig. 11. — Plots showing how an incorrect PA could affect the findings for the NGC 2366 data sets. The figure is formatted in the same 
way as Figure 9. From left to right are sets of plots for regions A, B, C, and T>, respectively. 



account for the difference shown for f2 and fl p in panels 
(h) and (i) of Figure 10. The solution method is clearly 
detecting something other than Q p ss Q. 

Figure 11 shows how an incorrect PA could affect the 
conclusions for the different regions of NGC 2366. For all 
four regions, the values of n do not change significantly, 
except for AMI at APA < -4°. Almost all of the values 
of D that are shown are too small to rule out the null hy- 
pothesis that the normalized residuals for AMI and AM2 
have a standard normal distribution, consistent with the 
values for these models in Table 3. The values of xt f° r 
regions A, B, and C, are fairly consistent, especially given 
their uncertainties. The values of F are all quite large 
for A and C. The smallest values of F occur for region 
V. Note that for all four regions, the values shown for 
the mean of fl p are within the 95% confidence bands for 
the value obtained using the adopted PA. It is reason- 
able to conclude that the accuracy of the adopted PA is 
not likely to significantly affect any of the conclusions for 
NGC 2366. 



5.3. Results for DDO 154 

Similar to NGC 3031, different regions: \y\ ^ 400", \y\ 
< 200", -400" < y < 0", and 0" < y < 400", produce 
consistent results. Therefore, only the results for most of 
the disk, \y\ ^ 400", are presented. This region provides 
24 independent calculations of Equation (1). 

The peaks in the distribution of Pks shown in Figure 
7 are not as narrow as they are in most of the other plots. 
Similar to most of the other plots, however, Pks is small 
near n = 0. Note that the value of Pks for the constant 
ftp model is the largest of the values of Pks shown in 
Table 3 for the constant £l p model. It is too large to 
rule out the null hypothesis that the distribution of the 
normalized residuals is a standard normal distribution 
for the constant tt p model. 

Plots of fL p calculated from the coefficients are shown 
in Figure 12. The material speed in the figure is calcu- 
lated from the rotation curve by Oh et al. (2011). All 
three models of tt p arc quite similar in the outer half of 
the disk. The plots for the constant fl p model and AM2 
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are consistent with each other to within their 95% confi- 
dence bands. The plot for AMI differs from the other two 
by showing corotation in the central part of the galaxy. 
None of the plots show corotation at ~ 100". 

The values of P F in Table 5 for AMI and AM2 are 
small enough to rule out the null hypothesis that a% = 0. 
This result for AM2 may be unreliable because the results 
of a Lllliefors test applied to its normalized residuals is 



inconclusive (Pl = 4.62%). Also note that the solution 
for this model has a value of xl that is smaller than 1 
by more than 2er. From the results for the F test, the 
results of the KS test for the constant fi p model, and the 
plots of Q p , the clumps of Hi in DDO 154 are corotating 
in the central and outer parts of the disk. Some of the 
clumps of Hi at ~ 100" are rotating more slowly than 
the material. The result for the constant il p model is the 
mean of Q p . 

Figure 13 shows how an incorrect PA could affect the 
conclusions for DDO 154. The results for DDO 154 are 
much more dependent on the value of the adopted PA 
than are the results for NGC 3031 or NGC 2366. The 
values of n,2 vary by as much as a factor of 6. The values 
of xt vary by almost as much, with the smallest values 
occurring around APA ~ 0°. There are values of APA 
where that the results of the F test are inconclusive, 
or that F is too small to rule out the null hypothesis 
that aj = 0. The results for the constant il p model 
deviates quite significantly beyond the 95% confidence 
intervals for the value obtained using the adopted PA. It 
must be concluded that the accuracy of the adopted PA 
affects some of the conclusions for DDO 154, particularly 
whether the solution using a constant Q p model, or one 
of the alternative models, provides a better solution for 
Equation (1). 

6. DISCUSSION 

6.1. How Different Properties of a Data Set Affect the 

Results 

The previous section shows the results from applying 
the solution method to data sets with different angu- 
lar resolutions, spectral resolutions, and sensitivities to 
structures on different scales. Results are also shown 
for galaxies without coherent or well-organized patterns. 
Knowing how these properties of a data set affect the re- 
sults is important for choosing data in future applications 
of the solution method. 

Of all the differences in the data sets, the angular res- 
olution had the most noticeable affect on the results by 
setting the number of independent calculations of Equa- 
tion (1). The most extreme cases can be seen in Figure 
8 for the NGC 3031 data sets. The results for the YP99 
data set demonstrate that a constant fl p , or the mean if 
there is evidence for a radial dependence, can be found 
with only 4 independent calculations of Equation (1), but 
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more calculations would be desirable in order to reduce 
the uncertainty. In order to determine whether there is 
a single, global pattern speed (i.e., Cl p is the same in dif- 
ferent regions of a galaxy), a number of independent cal- 
culations similar to what was used for AW96b and NGC 
2366 should be sufficient, given similar uncertainties. 

There is no noticeable affect on the results due to spec- 
tral resolution. Consistent results were obtained for the 
AW96a and H84 data sets despite a factor of ~ 4 dif- 
ference in the channel separation. Also consider that 
results with small uncertainties were obtained in Paper 
I for NGC 1365 using maps made from a 20.84 km s _1 
channel separation. The large channel separation for the 
NGC 1365 data may have been compensated for by hav- 
ing 62 independent calculations of Equation (1) for the 
whole disk. 

The results for the three different data sets of NGC 
3031 are remarkably consistent despite their different 
sensitivities to structures on different scales. The results 
for the AW96a and H84 data sets are almost identical 
even though they are more sensitive to small and large 
scales, respectively, and the AW96a data set is missing 
some flux. The results for the data sets of NGC 3031 
with the largest angular resolution, AW96b and YP99, 
are slightly less than the results for AW96a and H84. The 
smaller number of independent calculations of Equation 
(1) for the AW96b and YP99 data sets, and consequently 
larger uncertainties, may be the cause for this small dif- 
ference. 

The results for NGC 2366 and DDO 154 demonstrate 
that the solution method can produce meaningful results 
that are simple to interpretin the absence of a coherent 
or well-organized pattern. This shows that the solution 
method is not limited to grand-design patterns. It could 
be used to find il p for multiple armed and flocculent spi- 
ral galaxies as well. For these types of spiral galaxies 
the solution method would be measuring il p for different 
arms, or clumps of higher density material. Although 
this may violate the assumption of a well-defined pattern 
speed by TW84 (i.e., the surface brightness is assumed to 
be constant in a frame rotating at an angular speed O p ), 
the TW84 equations are generalizable to a frame that 
is rotating at a variable angular speed Ct p (r) (Engstrdm 
1994). Whether the results for galaxies without a coher- 
ent or well-organized pattern have any physical meaning 
can be determined using a variety of statistical tools such 
as xl an d the Lilliefor's test used in this paper. A very 
large or small value of xt > or a small enough value of Pl 
for the normalized residuals, could mean that Equation 
(1) is a poor description of the dynamics of a galaxy, or 
part of a galaxy, for example. 

This paper shows that having a sufficient number of 
independent calculations is important for applying the 
solution method. This limits the use of Hi radio synthe- 
sis data to nearby galaxies. For example, the resolution 
of the C-array of the VLA for Hi observing, which has a 
typical FWHM of ~ 13" for the synthesized beam, would 
be adequate for many nearby spiral galaxies. Radio syn- 
thesis data for tracers of molecular gas, which have a 
higher frequency, and thus a smaller synthesized beam 
(typically as small as sub arcsecond to a few arcseconds 
at ~ 100 GHz), would be useful for applying the method 
to more distant galaxies that have a smaller angular size. 



The resolution of optical data, such as Fabry-Perot Ha 
data, is typically ~ 1" - 2" due to atmospheric seeing. 
Such data would also expand the applicability of the so- 
lution method to galaxies with a smaller angular size. 

6.2. A Closer Look at the Findings for NGC 3031 

There is convincing evidence in the form of the F test 
that the spiral pattern in NGC 3031 is shearing. In fact, 
the values of F are quite large. This confirms the discov- 
ery of shear in the pattern by Westpfahl (1998). The pat- 
tern speed is very similar to the material speed, and there 
are no clear indications of unique corotation or Lindblad 
resonances. The same conclusions were found for the spi- 
ral pattern of NGC 1365. Note that the results of the 
F test are complimented by the results of the KS test 
for the constant f2 p model, from which it can be inferred 
that the constant £l p model is not the true model. It is 
important to point out that these findings for both NGC 
3031 and NGC 1365 are based on calculations that are 
independent of the assumptions made in spiral structure 
theories. 

The results for NGC 3031 are similar to the sheared 
gravitational instabilities of Goldreich & Lyndcn-Bell 
(1965) and the swing amplifier of Toomre (1981). Close 
passage of the companion galaxies NGC 3034 and NGC 
3077 could have caused the swing amplification in NGC 
3031. Several models agree that these galaxies passed to 
within 20 kpc of NGC 3031 about 600 Myr ago (Van der 
Hulst 1977; Cottrell 1976, 1977; Gottesman & Weliachew 
1975, 1977; and Brouillet et al. 1991). Simulations show 
that such interactions can produce grand design patterns 
that shear over time (S. Oh 2008; Dobbs et al. 2010; 
Struck et al. 2011). This implies that the spiral pat- 
tern in NGC 3031 is a transient feature that is driven 
by interactions with its companions, and in this way the 
winding dilemma is avoided. 

The results in this paper and in Paper I are part of a 
growing body of evidence that some grand-design spiral 
patterns are shearing, transient features. Some of this 
evidence is coming from simulations, which have yet to 
produce persisting patterns (e.g., Berman et al. 1978; 
Sellwood & Carlberg 1984; Sorensen 1985; Bottema 2003; 
Dobbs & Bonnell 2008; S. Oh et al. 2008; Dobbs et al. 
2010; Struck et al. 2011; Wada et al. 2011; Sellwood 2011 
for a review). Of the more recent simulations, Wada et 
al. (2011) snowed that the shearing, short lived (<~ 100 
Myr), recurring spiral arms in their simulations of a live 
disk have a pattern speed that is similar to the material 
speed. 

Simulations are useful tools for forming hypothesis 
about spiral structure, but a true understanding must 
ultimately come from observations. In Paper I it was 
pointed out that the original solution method of TW84 
produced results for at least 7 nearby spiral galaxies 
that are consistent with a shearing pattern. Foyle et 
al. (2011) also found evidence against long-lived, rigidly 
rotating patterns in spiral galaxies. They showed that 
the angular offsets between star formation tracers in 11 
nearby spiral galaxies, including NGC 3031, are not sys- 
temically ordered as predicted from a long-lived, rigidly 
rotating pattern. 

It is difficult to reconcile the findings in this paper for 
NGC 3031 with the density wave theory. All versions 
of the theory use as a foundation the premise that the 
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pattern is rotating approximately rigidly (Bertin & Lin 
1996, Chapter 2; Bertin 2000, Chapter 15; Binney & 
Tremaine 2008, Chapter 6). One possibility is that Q p is 
a step function, as advocated for by Meidt et al. (2008) 
and Meidt et al. (2009). The solution method used in 
this paper would not be able to detect such step functions 
because continuity is assumed for AMI and AM2. The 
forms used in Paper I, however, are capable of detecting 
such step functions, and did not find evidence for this 
form of Clp in the spiral pattern of NGC 1365. The same 
forms of ftp used in Paper I were also tried for the NGC 
3031 data sets, and there is no sign of a step function. 
Most importantly, note that there are no discontinuities 
in the spiral pattern of NGC 3031 as one would expect 
if Cl p is indeed a step function. 

7. SUMMARY 

In this paper the method derived in Paper I for solving 
the pattern speed equations of TW84 is applied to NGC 
3031, NGC 2366, and DDO 154 in order to show how 
different properties of a data set affect the results. The 
conclusions are as follows: 

1. ) Four different data sets of NGC 3031 produce con- 

sistent results despite differences in angular resolu- 
tion, spectral resolution, and sensitivities to struc- 
tures on different scales. They all produce convinc- 
ing evidence for shear in the pattern. The pattern 
speed is more similar to the material speed than it 
is to the speed of a rigidly rotating pattern. There 
are no unique points of corotation or 2-arm Lind- 
blad resonances. 

2. ) The angular resolution had the most noticeable 

affect on the results by setting the number of inde- 
pendent calculations of Equation (1). This is im- 
portant to consider when selecting data for future 



applications of the solution method, especially for 
determining the radial behavior of f2 p , and whether 
there is a single, global pattern speed. Lower angu- 
lar resolution data sets give few independent rows 
of data, so are less useful for determining Q p as a 
function of radius. 

3. ) The results for NGC 2366 show that the solution 

method can detect different pattern speeds in dif- 
ferent parts of a galaxy. 

4. ) The results for NGC 2366 and DDO 154 show that 

meaningful results that are simple to interpret can 
be obtained even in the absence of a coherent or 
well organized pattern. The results for these galax- 
ies also demonstrate that the solution method is 
not merely measuring the material speed, and thus 
serve as a counterexample to the results for NGC 
3031. 
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